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We use cosmic microwave background (CMB) and large-scale structure data to constrain cosmo- 
logical models where the primordial perturbations have both an adiabatic and a cold dark matter 
(CDM) isocurvature component. We allow for a possible correlation between the adiabatic and 
isocurvature modes, and for different spectral indices for the power in each mode and for their 
correlation. We do a likelihood analysis with 11 independent parameters. We discuss the effect 
of choosing the pivot scale for the definition of amplitude parameters. The upper limit for the 
isocurvature fraction is 18% around a pivot scale k = O.OlMpc -1 . For smaller pivot wavenumbers 
the limit stays about the same. For larger pivot wavenumbers, very large values of the isocurvature 
spectral index are favored, which makes the analysis problematic, but larger isocurvature fractions 
seem to be allowed. For large isocurvature spectral indices n- lso > 2 a positive correlation between 
the adiabatic and isocurvature mode is favored, and for rii so < 2 a negative correlation is favored. 
The upper limit to the nonadiabatic contribution to the CMB temperature variance is 7.5%. Of 
the standard cosmological parameters, determination of the CDM density ui c and the sound hori- 
zon angle 9 (or the Hubble constant Ho) are affected most by a possible presence of a correlated 
isocurvature contribution. The baryon density uib nearly retains its "adiabatic value". 



PACS numbers: PACS numbers: 98.70.Vc, 98.80.Cq 
I. INTRODUCTION 

The major part of the present cosmological data, 
including the cosmic microwave background (CMB) 
anisotropy and the large scale distribution of galaxies 
(large scale structure, LSS) is fit reasonably well by a 
simple cosmological model. This model has a spatially 
flat (fi = 1) background geometry. It has five energy 
density components, "baryons", photons, massless neu- 
trinos, cold dark matter (CDM), and a constant vac- 
uum energy (cosmological constant). The primordial 
scalar perturbations are gaussian, adiabatic, and scale- 
free (n a( j = const.). 

We call this model the "adiabatic model" in this pa- 
per. It has 5 parameters to be determined from the data, 
the Hubble constant, Hq = hlOQ km/s/Mpc, two den- 
sity parameters, cut = fli,h 2 and uj c = Q c h 2 (for baryons 
and CDM), and the amplitude A and spectral index n a( j 
of the primordial scalar perturbations. There is no ev- 
idence in the cosmological data for the presence of ad- 
ditional features or ingredients beyond this model, like 
tensor perturbations or neutrino masses, indicating that 
they are probably so small as not to show up in the data. 
(Actually there is also no evidence for a deviation from 
scale- invariance, n a( j = 1.) The concordance values of the 
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parameters are h ~ 0.7, ojb ~ 0.023, uj c ~ 0.12, n a( j ~ 1.0, 
and i~5x 10~ 5 . 

Besides these 5 fundamental cosmological parameters, 
there are 2 additional parameters needed when the mod- 
els are compared to CMB and LSS data: the optical 
depth t due to reionization, and the bias parameter b 
relating the observed galaxy power spectrum to the un- 
derlying matter power spectrum. 

The origin of the primordial perturbations is not 
known. The favorite candidate for their generation is 
quantum fluctuations during a period of inflation in the 
very early universe. While single-field inflation produces 
adiabatic perturbations, inflation with more than one 
field produces also entropy perturbations S in addition 
to the usual curvature perturbation 1Z. 

A general perturbation can be divided into an adia- 
batic mode + an isocurvature mode, where the adiabatic 
mode has no initial entropy perturbation, and the isocur- 
vature mode has no initial curvature perturbation. Al- 
lowing for the presence of an isocurvature mode does not 
improve the fit to the existing data (to the extent of jus- 
tifying the additional parameters), and thus there is so 
far no evidence for the existence of primordial isocurva- 
ture perturbations. However, it is of interest to find out 
what limits the data set to these perturbations, as the 
nature of primordial perturbations is an important clue 
to their origin. Moreover, the presence of an undetected 
isocurvature contribution may affect the determination 
of the main cosmological parameters. 

In principle, there can be different kinds of entropy per- 
turbations, and thus several different isocurvature modes. 
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Four different isocurvature modes were identified in 
the CDM and baryon isocurvature modes, and two neu- 
trino isocurvature modes. Allowing for the simultaneous 
presence of all four kinds would lead to so many parame- 
ters that it would be difficult to obtain meaningful results 
■ The signature of a baryon isocurvature mode in the 
data is rather similar to the CDM isocurvature mode, 
but weaker due to the smaller baryon density parameter. 

Here we consider only the CDM isocurvature mode 
in addition to the adiabatic mode. We allow the CDM 
entropy perturbations to have a different spectral index 
from the curvature perturbations, and to be (or not to 
be) correlated with them. In comparison to the adiabatic 
model this brings in 4 new parameters related to the am- 
plitudes and spectral indices of the entropy perturbations 
and their correlation with the curvature perturbations. 
Thus we have in total 11 parameters in our cosmologi- 
cal model. For sampling this 11-dimensional parameter 
space we use the Markov Chain Monte Carlo (MCMC) 
method. This is a follow-up paper of Q where a prelim- 
inary analysis (around the best-fit adiabatic model) was 
presented. Here we include more data in the analysis and 
sample the likelihood function more accurately. 

Before the Wilkinson Microwave Anisotropy Probe 
(WMAP) data became available, limits to the isocur- 
vature contribution in uncorrelated models had been ob- 
tained for the case n a d = fiiso = 1 in|5( and with n ac i and 
riiso as independent parameters in |(|, and in correlated 
models for one independent spectral index in Q. Pure 
CDM isocurvature models had been ruled out also in the 
case of a non-flat background geometry in p| . Correlated 
models were also studied in 

After WMAP, limits to correlated models were first 
obtained for the case of two independent spectral indices 
[ill H^ . In our earlier work Q we obtained prelimi- 
nary results for the case of three independent spectral 
indices using WMAP data only. Parkinson et al. 01 
considered a particular inflation model producing corre- 
lated perturbations. Moodley et al. [l4( considered mod- 
els with up to three isocurvature modes (CDM and two 
neutrino modes) present simultaneously, but all sharing 
the same spectral index. Ferrer et al. [l5j studied corre- 
lated perturbations resulting from inflaton and curvaton 
decay. They had two independent spectral indices. 

The most similar to the present study is that of Beltran 
et al. |16[. who consider one isocurvature mode at a time, 
and allow separate spectral indices for the adiabatic and 
isocurvature modes and their correlation. We compare 
their approach to ours at the end of Sec. [H] and their 
results to ours in Sec. IVIII 

Since the adiabatic and isocurvature components and 
their correlation are allowed to have different spectral 
indices, their relative amplitudes vary as a function of 
scale fc. We define the amplitude parameters at some 
chosen pivot scale fco- 

When the isocurvature component or the correlation 
is negligibly small, the corresponding spectral indices are 
not constrained by the data. Such conditionally uncon- 



strained parameters cause problems also for determining 
other parameters from the data. 

The way the isocurvature perturbation and correlation 
is parameterized (e.g. the choice of pivot scale fco) affects 
the integration measure in the parameter space. Thus 
different parameterizations correspond to different pri- 
ors. When parameters are weakly constrained by the 
data this ends up in different posterior likelihoods: When 
one parametrization (A) is used to obtain the likelihood 
function in the parameter space and the results are then 
expressed in another parametrization (B), the likelihood 
function is different from the case when parametrization 
B was used initially. (This difference can be "fixed" by 
importance weighting using the Jacobian of the param- 
eter transformation; but this does not address the ques- 
tion which parametrization is "correct" .) Such effects are 
discussed in Sec. IVII 

We find that the pivot scale should be chosen to be 
near the middle of the data sets used (in terms of lnfc). 
When the isocurvature spectral index rii SO is a free pa- 
rameter a wrong choice would spoil the analysis. This 
comes because the data does not prefer an isocurvature 
contribution. Then using fco that is close to the small fc 
end of the data, fc m i n , leads to extremely small (negative) 
n- lao , in order to minimize the isocurvature contribution. 
On the other hand, if fco is too close to fc max , then ar- 
bitrarily large isocurvature spectral indices are favored 
to minimize an overall isocurvature contribution in the 
range [fc m i n , fcmax]- Unfortunately, the "standard pivot 
scale" fc = 0.05MPC- 1 (used e.g. by CAMB [HEg) is 
quite close to fc max and another common choice (see e.g. 
[12)) to give for fco a value that corresponds to the present 
Hubble radius is nearly equal to setting fco = fc m in- 

When we started MCMC runs for our model we took 
fco = O.OSMpc" 1 , but soon realized that our Markov 
Chains ran towards artificially large rii ao . After fixing this 
problem, when we were finalizing the analysis of better 
runs with fco = O.OlMpc -1 , paper |16| with pivot scale 
fco = 0.05Mpc _1 came out. However, they had an ad hoc 
constraint n- lao < 3 that saved their main results from 
most of the artifacts that arise when the chains run to 
very large ni SO , With our choice of the pivot scale the 
likelihood for n- lso peaks at 7ij SO ~ 3 and drops rapidly 
around ni so ~ 4. Hence, the prior nj SO < 3 in 16] allows 
a comparison to our results. 

We obtain tight constraints to the CDM isocurvature 
contribution and find that, of the main cosmological pa- 
rameters, only the determination of lo c and h is sig- 
nificantly affected. Compared to the adiabatic models 
smaller values of u> c and larger values of h become accept- 
able when allowing for CDM isocurvature. Interestingly, 
although we have two additional degrees of freedom in 
spectral indices, dertermination of the baryon density is 
much less affected than in the models where all modes 
share the same spectral index. 

In Sec.^we introduce and motivate our parametriza- 
tion of correlated curvature and entropy perturbations. 
In Sec. IIHI we write down some technical details of our 
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MCMC study to determine these parameters, and in 
Sec. El we give and discuss our results. In Sec. \V\ we 
discuss the non-adiabatic contribution to the observed 
CMB and matter power spectra, and in Sec. I VII the ef- 
fect of changing the pivot scale. In Sec. IVIll we compare 
our results to those of [lfij . 



II. CORRELATED PERTURBATIONS 

The calculation of the CMB angular power spectra Ci 
and the matter power spectra P(k) starts from "initial" 
values lZ ra ,d and S ra d specified deep in the radiation dom- 
inated era (rad) , when all scales of interest are well "out- 
side the horizon" (i.e., the Hubble scale H -1 ). However, 
this "initial" time is well after inflation, or whatever gen- 
erated the perturbations, and refers to a time during and 
after which the evolution of the universe is assumed to 
be known. We denote the time when the perturbations 
were generated by the subscript *. For inflation, this cor- 
responds to the time when the scale in question "exited 
the horizon" (thus it is different for different scales k). 
Between t*(k) and i ra d the perturbation is outside the 
horizon, i.e., k is "superhorizon" . 

In the absence of entropy perturbations, curvature per- 
turbations remain constant at superhorizon scales. This 
is not, in general, true for entropy perturbations, which 
may evolve at superhorizon scales. Entropy perturba- 
tions may also seed curvature perturbations. This hap- 
pens, e.g., in two-field inflation, when the background 
trajectory in field space is curved [lfi Efl EH ] . 

Thus the relation between the "generated" and "ini- 
tial" values for TZ and S can be represented as Q 



n iad \ (i T ns \ (n t 

<S rad ) 10 T ss \S, 



(1) 



The transfer functions T xy (k) describe how the pertur- 
bations evolve from the time of inflation to the beginning 
of the radiation dominated era. The exact form of these 
functions is model dependent and that aspect is not stud- 
ied in this work. We approximate them by power laws. 

In the literature there are different sign conventions 
for the perturbations TZ and S. Wc define them so that 
an initial positive comoving curvature perturbation TZ la d 
corresponds to an initial overdensity 6 — Sp/ p > 0, and 
an initial positive entropy perturbation <S ra d corresponds 
to an initial CDM overdensity. In terms of the Bardeen 
potentials, $ and 4", defined so that the metric in the 
conformal-Newtonian gauge is 



ds 2 



(1 + 2$)dt 2 + a(t) 2 {l - 2^)6 ij dx i dx j , (2) 



the comoving curvature perturbation reads 
2 P ( 1 ^ , 5 



TZ = - 

3(p + p) \H dt 
and the entropy perturbation is 

S = S c - ^< 



(3) 



(4) 



where S c and S 7 are the CDM and photon density per- 
turbations. With these sign conventions, the ordinary 
Sachs- Wolfe effect is 



ST ~ 1 CP 

~ — g('M-ad 



2/ c 5 ra d) 



(5) 



(here f c = Lu c /(^b + w c )), and a positive correlation be- 
tween 7?. ra d and <S ra d leads to an additional positive con- 
tribution to the large scale CMB anisotropy, and also to 
a positive contribution to the matter power spectrum. 

We define the correlation C xy (k) between two pertur- 
bation quantities (random variables), x and y, as 



(x(k)y*(k'))\ tad = ^C xy (k)6^(k-k'). 



(6) 



The transfer function T-jis(k) leads to a correlation 
between TZ ra d and <S ra d from uncorrelated TZ* and <S*, 

Cnn{k) = V n (*, k) + T ns {k) 2 T s {*, k) (7) 
Cns(k) = T ns (k)T ss {k)V s (*, k) (8) 
C ss (k)=T ss (k) 2 V s (*,k), (9) 

where Vn{*, k) and Vs(*, k) are the power spectra of TZ* 
and <S*. 

Approximating the power spectra V-jz(*,k), Vs(*,k) 
and the transfer functions T-jis(k), T$s(k) by power laws 
with spectral indices mi, mi, m§, and 7714, respectively, 
we get that the autocorrelations (power spectra) have the 
form 



T n {k) ~Cnn{k) = A 2 , 



V s {k)=C ss {k) = B 2 [- 



"adl-l 



■At 



ko 



»ad2 — 1 



(10) 



where n a di = mi + 1, n a d2 = ^2 + 2m3 + 1, and 
ni S o = TO2 + 2mA + 1 and the epoch (rad) is implied. The 
three components are the usual adiabatic mode, a second 
adiabatic mode generated by the entropy perturbation, 
and the usual isocurvature mode, with amplitudes A r , 
A s and B at the pivot scale k = ko, respectively. 

The cross-correlation between the adiabatic and the 
isocurvature component is now 



Ciis{k) — CsTz(k) — A S B 



(11) 



where n COI = mi+m^+mi+l = (rii S o+«ad2)/2. The cor- 
relation is between the second adiabatic and the isocur- 
vature component as is natural since these components 
have the same source. 

We have chosen the pivot scale ko — 0.01 Mpc" 1 , 
but we also consider pivot scales 0.002 Mpc -1 and 
0.05 Mpc -1 in Sec. I VII We shorten the notation by defin- 
ing k = k/ko- 
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The CMB angular power spectrum is given by 



(T/E/B)(T/E/B) 



4 -E / f C ^T /B) (k) 9 ^ B \k), (12) 

xy J 

where the gi's are the transfer functions that describe 
how an initial perturbation evolves to a temperature (T) 
or polarization (E- or B-mode) anisotropy multipole 

Now, using the equations (flTl)) . (fTTJl and l(T2"|l we obtain 
for the temperature angular power spectrum 



2A s Bgl l gl l k^- 1 



2/ s »TTadl 



--A r C l 



2/YTTad2 



^2^iTTiso 



A s BCr mT - 



(13) 



and for the TE cross-correlation spectrum 



C, TE - 47T 



dk r 

i D 2 T E fniso-1 

+ -B 9si9si k 

+ A s B(gl l gV l +gJ l gl l )k n ^- 1 

_ / |2/.TEadl i / i2/.TEad2 , D 2/»TE 



^r5^ zffTC ;^" adl 1 + ^9tz i9tz ;^™ a 



A<BC? Ec 



(14) 



There are thus 3 amplitude parameters (three abso- 
lute values and one sign, the relative sign of A s and B.) 
Now we need to choose the amplitude paramctrization 
to be used in the likelihood analysis, i.e., what shall we 
use as the three independent parameters with flat prior 
likelihoods. One choice would be just A r , A s , and B. 
However, we would rather express our results in terms 
of a total amplitude, a relative isocurvature contribution 
and a correlation. 

In Q we followed an d used 



/iso — 



B 2 



A 2 -j- A 2 

for the isocurvature contribution and 



e [0, oo) 



cos A ee sign(A s _B) 1 



A 2 , 



(15) 



(16) 



for the correlation (with the sign convention opposite to 
that of 

EH)- 

The data is quadratic in these parameters 
(see Eq. (E3J), meaning that fairly large values of /i so 
and cos A are needed for the effect to show up in the 
data. This exacerbates the problem that models with 
a small A s and B get a lot of weight in the likelihood 
function, since the spectral indices n a( j2 and n iso are not 
constrained. 



A flat prior for cos A leads to a non-flat prior distri- 
bution for cos 2 A. Thus the parametrization by cos A 
favors small multiplier cos 2 A in front of the second adi- 
abatic component in 0, [n|. Moreover, large values of 
sin 2 A = 1 — cos 2 A are then favored, so that even without 
any data the first adiabatic component will be favored in 
the likelihood analysis. Likewise, the parametrization by 
/i so (instead of f 2 so ) favor small multiplier in front of the 
isocurvature component. All in all, there was an implicit 
bias towards pure adiabatic models in HE3. A similar 
caveat applies to [l5j . 

We would prefer amplitude parameters for which the 
data has a linear response. We define a total amplitude 
parameter A by 



A 



Ai + A 2 



B 



(17) 



and the isocurvature fraction and correlation parameters 

(18) 



B 2 

a ee — 6 [0, 1] 



7 ee sign(A s B)- 



A 2 



e[-i,l] 



(19) 



A 2 + A 2 

Now the total angular power spectrum can be written as 
C, = A 2 [(l a)(l - | 7 |)CT dl + (1 - a)| 7 |Q ad2 



+ aCf + sign( 7 )Va(l-a)| 7 |C*r] 



(20) 



Here Cf dl and Cf d2 represent adiabatic spectra which 
would result from a curvature perturbation 7?. ra d with 
unit amplitude (A r = 1 or A s = 1) at the pivot scale kg. 
(They are otherwise the same, but have spectral indices 
n a( ji and n a d2-) Likewise, Cf° represents an isocurvature 
spectrum from a CDM entropy perturbation of unit am- 
plitude (B — 1), and Of 01 the extra contribution from 
correlation for A S B = 1. (See Figs, ^and^ which rep- 
resent the case of scale- invariant perturbations.). The 
"hatless" C z adl , Cf d2 , Cf°, which are necessarily non- 
negative, and Cf or , which can also be negative, are the 
contributions to the total C/. A relation similar to l|2(J|) 
holds for the matter power spectrum P{k). 

Note that, e.g., a = 0.5 does not mean that the adi- 
abatic and isocurvature contributions would be equal at 
any particular scale. Since a refers to the ratio of pri- 
mordial perturbations, to which the C/ contributions are 
related through the transfer functions, the situation is 
different for different scales, and depends on the other 
cosmological parameters. In particular, if the spectral in- 
dices are very different, a very small isocurvature fraction 
can still correspond to a large isocurvature contribution 
at some scales and vice versa. 

Wc define a shorthand notation 



a CO r = si; 



;n(7)Va(l - a)\j\ 



(21) 



for the relative "weight" of the correlation spectrum 
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Figure 1: The unit-amplitude component angular power spec- 
tra Cf d (red), Cf° (blue), and C^ 1 (green) of Eqs. O and 
12Ut for the case of spectral indices n a d = ni so = 1 and other 
cosmological parameters representing median values of their 
marginalized likelihoods from our 11-parameter model. These 
curves would represent the relative contributions to the total 
Ci for the case a = 0.5, 7 = 1, i.e., "equal" weights for the 
adiabatic and isocurvature contributions and a maximal pos- 
itive correlation between them. 



The problem remains that when some multiplier in 
(|20|l is close to zero, the spectral index of the correspond- 
ing component becomes unconstrained leading to more 
volume in parameter space upon marginalization. This 
may introduce a bias towards "pure" models where the 
isocurvature or correlation amplitude is zero. 

We want the pivot scale to be roughly in the middle of 
the data set used, and have chosen fc = O.OlMpc" 1 as 
our pivot wavenumber. For the concordance values of the 
cosmological parameters, fl\ = h = 0.7, this corresponds 
to "pivot multipole" Iq ~ 140. [The correspondence is 
Iq ~ D^ko, where = .D*(/i, Q m ) is the angular 
diameter distance to last scattering. (h, 0.7, 0.3) « 
h^lO OOOMpc -1 while the "old day's standard value" 
was D*(h, 0, 1) w h~H OOOMpc" 1 .] 

This work is similar to a recently published study by 
Beltran et al. [16|, The main differences are: 1) Dif- 
ferent parametrization of correlation. When we divide 
the adiabatic spectrum in a correlated and an uncorre- 
cted part, they consider the total adiabatic spectrum Vn 
and the correlation spectrum C-rs as the basic entities, 
which they approximate by power laws. This leads to 
constraints on the correlation spectral index n cor , which 
depend on the correlation amplitude, and therefore they 
introduce a related parameter, ll S cor " , to be the indepen- 
dent parameter, leaving derived parameter. 2) 
They have set an upper limit ni SO < 3, whereas we allow 
niso to vary over a wider range. 3) They use a pivot scale 
k = 0.05 Mpc" 1 (l - 700). We use k = 0.01 Mpc" 1 
(Iq ~ 140), but consider also the effect of changing the 
pivot scale. 4) They use a larg er data set, including type 
la Supernova (SNIa) data [22] , whereas we use CMB and 
LSS data only. 5) They include an equation-of-state pa- 




Figure 2: The same as Fig. JTJ, but for (a) (7, TE and (b) the 
matter power spectrum P(k). We also show (c) the (7 ; TT of 
Fig. Q with a logarithmic scale, so that the effect of chang- 
ing the spectral indices can be readily estimated from the 
figure. The pivot scale fco = O.OlMpc -1 becomes ko/h — 
0.01418Mpc _1 for the parameter values used (h = 0.7053) for 
this plot. 



rameter w for dark energy, wheras we keep w = — 1. 6) 
They consider neutrino isocurvature modes also. 

Crotty et al. 01 an d Beltran et al. Hi use the same 
isocurvature parameter a as we use, but they use the 
correlation parameter 



cos A =— sign(A s B)i 



-signWVM 
(22) 
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In ^3 P is assumed scale invariant, whereas in 0] it is 
approximated by a power law with index 

f^QQT SO that our 

= 2( n ad2 +^iso) corresponds to their n c 

^iso) ■ 
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III. TECHNICAL DETAILS OF THE ANALYSIS 

The model we are studying has 11 parameters. We 
have chosen to use the following independent parame- 
ters (primary parameters) for the likelihood analysis: the 
baryon density u>b, the CDM density u c , the sound hori- 
zon angle 8, the optical depth due to reionization t, the 
bias parameter b, the uncorrelated adiabatic spectral in- 
dex Uadij the correlated adiabatic spectral index n a d2, 
the isocurvature spectral index nj SO , the logarithm of the 
overall amplitude ln(10 10 A 2 ), the isocurvature fraction a 
and the correlated fraction 7 of the adiabatic perturba- 
tions. 

The sound horizon angle (in units of ^55 radian) 



6(u b ,uj c ,h) = 100 x 



/' 



(23) 



where s* is the sound horizon at last scattering and -D* 
is the angular diameter distance to last scattering 23] , is 
used as an independent parameter instead of h (or f2/v), 
since it is more tightly constrained by the data. 
The bias b is defined by 



jSDSS 



Zeff^0.15 



= b 2 P(k) 



(24) 



So we multiply the present-day theoretical matter power 
P(k) by b 2 before comparing to the galaxy power spec- 
trum observed by the Sloan Digital Sky Survey (SDSS) 
[24| at effective redshift z e g. In the figures, we actually 
plot b 2 P(k). 

We find the posterior likelihoods for the primary pa- 
rameters and a number of derived parameters using the 
Markov Chain Monte Carlo (MCMC) method. The 
chains are generated using our modified version of the 
publicly available CosmoMC code [13 • The CMB an- 
gular power spectra and the matter power spectra are 
calculated by the CAMB code [II S3 (see also (H). It 
needed some modifications for faster treatment of corre- 
lation. 

CosmoMC/CAMB evaluates the matter power spec- 
trum in linear perturbation theory. However, very small 
scales (k/h > 0.15Mpc _1 ) have already become non- 
linear. The publicly available code HALOfit utilizes re- 
sults from lattice simulations of clustering [27|. How- 
ever, the applicability of the HALOfit to our model is 
not granted, since the lattice simulations have been per- 
formed in adiabatic models with moderate spectral in- 
dices. Hence, following the recipe of [24j, we calculate 
the matter power spectra in linear theory and compare 
them only to the first 17 data points (k/h < 0.15Mpc _1 ) 
of the SDSS galaxy survey [2|. 



For the observational CMB data we take the WMAP 
temperature autocorrelation (TT) and temperature- 
polarization cross-correlation (TE) data [23, 123, H3| • To 
extend the coverage of the data to higher multipoles we 
use the TT data from CBI [p and ACBAR (33, which 
we later call "other CMB data" . 

Details of the data sets are: 

• WMAP TT, 899 data points, 1 = 2 - 900, 
(k ~ 1.4 x lO^Mpc" 1 - 0.064Mpc _1 ). 

• WMAP TE, 449 data points I = 2 - 450, 
(k ~ 1.4 x lO^Mpc" 1 - 0.032Mpc _1 ). 



• ACBAR TT, 7 ^-bands, l cS = 991 
(k ~ 0.071MPC" 1 - 0.131MPC- 1 ). 



1831, 



• CBI TT, 13 Z-bands, l eS = 369 - 1884, 
(k ~ 0.026MPC- 1 - 0.135MPC- 1 ). 



• SDSS galaxy power, 
0.016MPC" 1 - 0.15Mpc 



17 fc-bands, k e g/h 



In parenthesis we indicate what wave numbers the given 
multipole ranges correspond in models that have fl\ = 
h = 0.7, i.e. D» ~ MOOOMpc" 1 . The total number 
of data points (1385) leads to the reduced number of 
degrees of freedom v = 1385 — 11 = 1374 for our model 
and v = 1385 — 7 = 1378 for the adiabatic model. 

First we did several 8-chain runs to see what happens 
in a MCMC study of our model. Finally, we chose a suit- 
able parametrization, described above, and performed an 
8-chain initialization run with the option to update the 
proposal matrix (jump function) turned on in CosmoMC. 
We used this run to obtain a good proposal matrix for our 
full run. In our full run we ran the code on an IBM AIX 
cluster utilizing 32 processors for 12 days to produce 32 
chains that started from separate randomly picked points 
in parameter space. After cutting off the burn-in periods 
the total number of accepted steps, i.e., different com- 
binations of our primary parameters, was 266 651. The 
total number of different models tried (step trials) was 
8 005 143. The option to update the proposal density 
while generating the chains was not used in order to pro- 
duce pure MCMC chains. In addition to this main run, 
another set of 8 chains with 60 254 different models with 
continuously updated proposal density is used as addi- 
tional data when discussing the effect of the pivot scale 
in Section IVII For a clear review of steps included in 
MCMC analysis, especially the meaning of marginalized 
likelihoods, see the Appendix of Tegmark et al. |33|. 

The parameters were allowed to vary within the fol- 
lowing ranges: 

cu b G [0.005,0.1], uj c e [0.01,0.99], 6 £ [0.3,10.0], 

n a diG[-3,4], n ad2 e[-3,4], n iso G [-3, 12], 
t e [0.01, 0.3], ln(10 10 A 2 ) S [1, 7], 66 [0.1, 2.5], 
7 G [-1, 1] and a G [0, 1] 

The MCMC method implicitly assigns flat priors for 
these independent parameters. The ranges for a and 
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7 follow from their definitions. For the other parame- 
ters, except r, we have set very wide ranges, so that the 
likelihood is negligible at the boundaries. However, we 
also imposed a top-hat prior for the Hubble constant: 
0.4 < ft < 1.0, which cuts off some models that would 
otherwise be acceptable (at 95 % C.L.). 

We have constrained r to be less than 0.3. We found 
in our preliminary studies that there are models with 
r > 0.3 that fit well to the data. These models form a 
separate region in the parameter space, and have also a 
high baryon density, of the order of o>b ~ 0.03. This high 
baryon density is much above the values obtained from 
big bang nucleosynthesis (BBN) calculations [3J( and we 
decided not to consider such models in this paper. In- 
cluding the t > 0.3 region would be problematic with the 
MCMC method as it is not well suited for such bimodal 
distributions. Moreover, r > 0.3 leads to a very high 
reionization redshift, which is not favored by astrophysi- 
cal considerations |35| . 

To cover our parameter space as well as possible, 
within the limits of available computational resources, 
the starting point for each of the 32 chains was randomly 
selected from the following Gaussian distributions: 





= 0.0236 ±0.005, 


L0 C 


= 0.124 ± 0.035, 





= 1.047 ± 0.038, 


T 


= 0.11 ±0.229, 


™adl 


= 0.97 ±0.27, 


«ad2 


= 0.97 ± 1.60, 


<Hsa 


= 2.09 ±3.70, 


b 


= 0.99 ±0.34, 


7 


= 0.01 ± 1.3, 


a 


= 0.035 ±0.24, 


bx(10 w A 2 ) 


= 3.20 ±0.4. 







The width for a given parameter is four times the width 
of the posterior distribution of the same parameter from 
our preliminary runs. 



We first discuss the effect of allowing a (possibly corre- 
lated) isocurvature contribution, on the determination of 
the standard cosmological parameters. The likelihoods of 
Uf,, uj c , 9, n a di, ln(10 10 A 2 ), t, and b, are compared with 
the corresponding likelihoods of the adiabatic model in 
Fig. El 

The amplitude A has now a different meaning than in 
the adiabatic model, as it includes the isocurvature con- 
tribution also. Since the isocurvature transfer functions 
lead to less power in most of the data from a given pri- 
mordial amplitude than the adiabatic transfer functions 
(see Figs.n an d|2J>, larger total amplitudes A are allowed 
for models with a significant isocurvature contribution. 

The distribution for the adiabatic spectral index n a( ji 
has become much wider. The reason for this is that the 
correlated adiabatic component ("ad2") may take the 
role of the adiabatic perturbation of the adiabatic model: 
If |7| ~ 1, but a is small, the model looks like the adia- 
batic model; that the adiabatic mode is correlated with 
the isocurvature mode does not have much significance, 
if the isocurvature component itself is negligible. In this 
case n a d2 is then constrained to be close to the spectral 
index value of the adiabatic model, but n a( ji becomes 
unconstrained, as this contribution has negligible ampli- 
tude. We discuss the question of the adiabatic spectral 
index further in Sec. IIV Al 

The uncertainties in the determination of lu^, t, and 
b are increased somewhat. We discuss lo c and 6 in 
Sec HVBl and u b in Sec. HVT1 We devote Sees. HVHl 
and IIV El to the isocurvature and correlation parameters, 
respectively. In Sec. HVFl we discuss a warning example 
of a model with very large n- mo that must be rejected for 
several reasons. 

A. Adiabatic spectral index 



IV. RESULTS 

In Fig. we show the marginalized ("1-d") likelihoods 
for those 7 of our independent parameters, which corre- 
spond to the 7 parameters of the adiabatic model. In 
Fig. 01 we show likelihoods for some derived parameters 
related to them. 

In Fig. [5] we show the marginalized likelihoods for our 
remaining 4 independent parameters, and in Fig. [S] for 
some related derived parameters. 

Flat priors for our independent parameters lead to non- 
flat priors for the derived parameters, which contribute 
to some features in the distributions of the latter. 

The best-fit (11-parameter) model has % 2 = 1459.29, 
just slightly better than the best-fit (7-parameter) adi- 
abatic model \ 2 — 1459.65. Thus there is clearly no 
indication in the data for the presence of an isocurvature 
contribution. Our results should be considered in terms 
of upper limits to isocurvature perturbations and uncer- 
tainties in the determination of cosmological parameters 
due to the possibility of an isocurvature contribution. 



We can define an effective single adiabatic spectral in- 
dex by 



dlnV n {k) 
dlnfc 



= Kdi - i)(i - MJfc"- 1 - 1 + Kd2 - lJMfc"-'- 1 

(1 - |7|)fc«adl-l + | 7 |fc" a d2-l 

which is scale dependent. The first derivative 

dnllCk) = (l-|7l)l7lKdi-n ad2 ) 2 fc"^+"^ 
d\nk [(1 - |7|)A: n -ii + | 7 |fc«ad 2 ]2 

is zero only when n a di = ra a d2 or 7 = 0, ±1. Otherwise 
it is positive 

At our pivot scale we have k — 1 and the above ex- 
pressions simplify to 

'ad| fc=fe0 - Kdi - 1)(1 - |7l) + (n ad 2 - l)|7l + 1 , (25) 



and 



dnim 



dhi k 



k—kn 



(«adl-«ad 2 ) 2 (l-|7l)l7|- (26) 
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Figure 3: Marginalized likelihood functions for the standard 
cosmological parameters (i.e. those that exist in the adiabatic 
model). The solid (black) line is the likelihood in our 11- 
parameter model, the dashed (red) line is for the adiabatic 
model. Other line types show the effects of additional priors 
discussed in the text: dotted (blue) for Gaussian Qa = 0.70 ± 
0.04, and dot-dashed (green) for Gaussian 7 = 0.0 ± 0.02. 



From Fig. |3| we observe that n ac ii is much more loosely 
constrained than the n a a of the adiabatic model. The dis- 
tribution for n a( j2 becomes even wider than the one for 
rt a di, see Fig. [S] The reason is that the MCMC chains 
contain many models with I7I close to zero allowing ro a d2 
to take any value or | | close to 1 allowing n a di to take 
any value. However, the effective adiabatic spectral index 
(I25|l becomes nearly as tightly constrained as the spectral 
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Figure 4: Marginalized likelihoods for two derived parame- 
ters, Qa and Ho- The line styles have the same meaning as 
in Fig. 01 
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Figure 5: Marginalized likelihoods for the parameters related 
to the isocurvature mode and correlation. We show the 4 
remaining independent parameters, a, 7, ni so , n a d2- The solid 
(black) line is the full likelihood, other line types show the 
effects of additional priors discussed in the text: dotted (red) 
for Gaussian Qa > 0.82, dashed (blue) for Gaussian Qa = 
0.70±0.04, and dot-dashed (green) for Gaussian 7 = 0.0±0.02. 



index in pure adiabatic models. The 95% C.L. regions 
are 0.910 < n|f < 1.050 with median 0.968 and 0.923 < 
n a d < 1.013 with median 0.961, see also Fig.[7fa). More- 
over, the data disfavor (positive) running of the adiabatic 
spectral index. For the 95% C.L. upper limit we obtain 
dn^/dlnk < 0.03 at k = 0.01Mpc~\ see Fig. ^b). The 
largest k in the data sets is about fc m ax ~ 0.15Mpc _1 . 
So the maximum running from fco to fc max is approxi- 
mately An = 0.03 x ln(fc max /fco) = 0.08. The quadrupole 
(/ = 2) corresponds to k m i n w 1.4 x 10 _4 Mpc _1 leading 



Figure 6: Marginalized likelihoods for the isocurvature- 
related derived parameters, a C or and n CO r. The line styles 
have the same meaning as in Fig. 
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Figure 7: (a) Marginalized likelihoods for the effective adia- 
batic spectral index (solid line, black) compared to the 
spectral index of the pure adiabatic model n a( j (dashed line, 
red), (b) Marginalized likelihood for dn^/dlnk at the pivot 
scale k = 0.01 Mpc" 1 . 



to An = 0.03 x ln(fc min /fco) = -0.12. 
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Figure 8: (a) 2-d marginalized likelihood for and lo c , show- 
ing that the large sound horizon angles 9 are connected with 
low CDM densities u c . We indicate the 68% (solid) and 95% 
(dashed) C.L. regions for our isocurvature model (black) and 
for the adiabatic model (red), (b) 2-d likelihood for the two 
derived parameters Ho and Q,a closely related to the indepen- 
dent parameters 8 and lo c . 



B. Small matter density models 

In Fig. |3 the most obvious difference from the adia- 
batic model is the extension of the 9 likelihood towards 
larger sound horizon angles and the uj c likelihood towards 
smaller densities. These two features are related as can 
be seen in Fig.^a). The corresponding effect is seen in 
the two derived parameters, f^A, Hq, closely related to 9 
and uj c , see Fig.|HIb). Compare to a simalar Fig. in |36j . 

The 1-d likelihood for the derived parameter show 
(Fig. a second peak at fl\ ~ 0.87. (This feature is 
somewhat enhanced because the flat prior for our inde- 
pendent parameters actually leads to an increasing prior 
for the derived parameter Q\, and larger values of fl\ 
are cut off with our h < 1 constraint.) 

Thus the possibility of an isocurvature contribution 
leads to larger J7a models becoming acceptable by the 
CMB and LSS data. According to Fig. ED these mod- 
els have a positive correlation between the adiabatic and 



isocurvature modes. Indeed, if we cut to the subset of 
"uncorrelated models" , 7 = 0.0 ±0.02, the (large 9, small 
u) c ) feature disappears from the 1-d likelihoods. 

In Figs. and HJ we show the 1-d likelihoods of the 
isocurvature-related parameters separately for the large- 
f^A subset (dotted red lines) and with a 7 = 0.0 ± 0.02 
prior that cuts the more correlated models off (dot-dashed 
green lines). We see clearly that the large- ft\ models 
are associated with a positive correlation between the 
isocurvature and adiabatic modes. 

The angular and matter power spectra of the best-fit 
large-ilA model (from the subset Qa > 0.82) are shown 
in Fig This model has x 2 = 1461.86. Compared 
to the best-fit adiabatic model, the somewhat worse fit, 
Ax 2 = 2.21, is due to 1) a worse fit to the SDSS data 
(Ax 2 = 1.83) and 2) a worse fit to the Sachs- Wolfe region 
(2 < I < 21) of the WMAP TT data (Ax 2 = 2.51). The 
latter is due to the increased late ISW effect caused by 
the larger This model fits the rest of the CMB data 
better than the adiabatic model. 
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Figure 9: 2-d marginalized likelihood for a C or and SIa, show- 
ing that the large SIa are connected with a positively corre- 
lated isocurvature contribution. 



The reason the larger sound horizon angles (which shift 
the acoustic peaks left, i.e., towards smaller I) are ac- 
cepted is the correlation contribution C ; cor , whose acous- 
tic peaks are at somewhat larger I than the adiabatic 
ones, and thus adding it appears as a shifting of the 
peaks to the right (i.e., towards larger I). An uncorre- 
cted isocurvature contribution cannot do the same trick, 
since the isocurvature acoustic peaks are too much to the 
right for adding them to appear as a shift in peak posi- 
tion. The distribution of the isocurvature spectral index 
is concentrated at the upper end of the allowed range for 
n iso in these S}\ > 0.82 models, see again Figs. |S] and 
This is required for the correlation contribution to 
maintain roughly the same relative power through the 
acoustic peak region. 

Because of this large spectral index, especially the cor- 
relation contribution also changes the shape of the matter 
power spectrum, see Fier. llOf b). This allows for a smaller 
"shape parameter" Sl m h to fit the SDSS data, than the 
SDSS result Q m h = 0.21 ± 0.03 for adiabatic models. In 
the adiabatic model, large values of SIa and h would be 
allowed by either the CMB or the LSS data alone, but 
not by the combined data sets, because either data set 
allows a narrow region (the "vanilla banana" in Fig. 5 of 
|3^ ) in the (Q\,h) plane, but these regions have some- 
what different orientations. The correlation contribution 
makes both regions wider, in such a way that their over- 
lap is extended to higher h and smaller Sl m (larger S}\), or 
in terms of our independent parameters, towards smaller 
u> c . In fact, even ft > 1.0, with Sl m < 0.1 (or Sl\ > 0.9) 
would be allowed, but our ft < 1.0 prior cuts them off. 
These models also favor smaller bias parameters b and 
baryon densities 

One might expect the above to work in the other di- 
rection too, negative correlation allowing models with a 
smaller 9, SIa, ft, and a larger w c , but apparently some 
other feature in the data prevents the larger u c required. 

The large values of f^A are ruled out by the high-z 
Type la Supernovae redshift-magnitude (SNIa) data|22j. 




k/h [Mpc 



Figure 10: (a) The TT spectrum for the best-fit large SIa 
(SI a > 0.82) model, showing how the correlation part con- 
tributes towards shifting the acoustic peaks to the right. This 
model has u b = 0.0221, lo c = 0.100, 9 = 1.082, r = 0.177, 
b = 0.796, ln(10 10 A 2 ) = 3.24, n adl = 0.949, a = 0.108, 
fiiso = 3.11, 7 = 0.57, n a d2 = 1.043. The correlated adiabatic 
component ( "ad2" , magenta) dominates over the uncorrelated 
adiabatic component ("adl", red). The other curves are the 
correlation component ("cor", green) and isocurvature com- 
ponent ("iso", blue). The total angular power ("tot", black) is 
a sum of these components, (b) The matter power spectrum, 
showing how the correlation contribution changes its shape. 



Therefore the large- SIa feature was not seen in ^(|. We 
did not use the SNIa data; but to study the effect of a 
SNIa constraint we simulated it by importance weighting 
our MCMC chains with a Gaussian fi A = 0.70 ± 0.04 
distribution. We show the 1-d likelihoods both with and 
without this extra prior in Figs. [3HHI The effect of this 
SNIa constraint cutting the (large 9, small u c ) models off 
is clearly seen in them. 



C. Baryon density and Hubble parameter 

In pure adiabatic models the baryon density is prac- 
tically determined by the heights of the first and sec- 
ond acoustic peaks (and the valley between them). An 
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isocurvature contribution modifies these heights and thus 
one expects looser constraint for uj b in mixed adiabatic 
and isocurvature models. However, our constraints 0] 
(0.0220 < oj b < 0.0246 at 68% C.L., 0.0207 < uj b < 
0.0263 at 95% C.L., median 0.0232) are very close to the 
adiabatic model (0.0221 < uj b < 0.0240 at 68% C.L., 
0.0213 < uj b < 0.0250 at 95% C.L., median 0.0230). 
Moreover, we have checked that the isocurvature ampli- 
tude (a) dependence of the constraints for us b is very weak 
within the allowed range a < 0.18. 

As can be seen in Fig. the median of uj b shifts only 
marginally towards larger values regardless of the (extra) 
priors chosen. Our result is consistent with 0] where 
the 1-d likelihoods for the adiabatic reference model 
and for the correlated CDM isocurvature model were 
practically indistinguishable. (The neutrino isocurvature 
modes shifted u> b towards smaller values unlike in some 
other models to be discussed below.) 

Both our result and the likelihood for ui b in 16] differ 
from ^4| where the median shifted significantly towards 
larger values (0.023 < oj b < 0.029 at 68% C.L., median 
0.026). In their model the adiabatic, CDM isocurvature 
and correlation spectral indices were ke pt e qual, i.e. they 
had only one free spectral index. Both |16| and used 
CMB and LSS data sets very similar to those used by 
us. On the other hand, in [15| the curvaton decay calcu- 
lation (see e.g. 01) was extended to the case when the 
curvaton does not necessarily behave like dust. The re- 
sulting correlated CDM isocurvature perturbations from 
the mixed inflaton-curvaton decay (or e.g. from double- 
inflation which produce primordial power spectrum of 
similar type) were considered in the light of WMAP data 
alone. Then uj b got even larger values. The 68% C.L. 
region obtained in 15j] was 0.027 < uj b < 0.042 with 
median 0.032. (The best-fit model had uj b = 0.041.) Al- 
though, the form of the primordial power spectrum in 
[TEllsiq l looks quite similar to ours, the important differ- 
ence is that the adiabatic and isocurvature components 
have equal spectral indices there. Note however, that as 
the curvature and entropy perturbations are both a sum 
of two components the model has two independent spec- 
tral indices insted of just one of 0. In [3^ a model 
with equal spectral indices for adiabatic, CDM isocurva- 
ture and neutrino isocurvature modes yielded also very 
large ui b ~ 0.04. The CMB data alone led to a bit larger 
w b than CMB and LSS data together. 

Three years ago Trotta, Riazuelo, and Durrer demon- 
strated in ^3 that allowing for "general isocurvature 
modes" (adiabatic, CDM isocurvature and neutrino 
isocurvature with equal spectral indices in their study) 
prevented one from obtaining an upper bound for u b from 
that day's CMB data (COBE HHlUm and Boomerang 
p4|). In ^(j most of the new freedom for u b was ex- 
plained to come from the neutrino isocurvature density 
mode which can adjust the height of the second acoustic 
peak more than other isocurvature modes. Hence, one 
would expect more freedom for Lu b when allowing also 
for a neutrino isocurvature density mode instead of just 



a CDM isocurvature mode (or a neutrino isocurvature 
velocity mode). However, one can not see this effect in 
(l(| where the median of iv b was shifted a bit towards 
smaller values (compared to other cases) and the width 
of the 1-d distribution remained small. While part of the 
different effect of the neutrino isocurvature density mode 
in [4(j and [0| could result from the different data sets 
used (the former used CMB only, the latter used pre- 
cision CMB and LSS), we think that the fundamental 
explanation resides in spectral indices. The same applies 
to other isocurvature modes. 

To be explicit, our model and the model in 

[H (both 

have three independent spectral indices) yield very lit- 
tle difference to the adiabatic case whereas models stud- 
ied in [H HI El (all have "n ad = n iso ") lead to 
larger medians and wider distributions for uj b . While 
we have not studied in detail the reason for this differ- 
ence, we discuss one possibility. The CMB data forces 
the dominant adiabatic spectrum close to scale invari- 
ance (n a( j ~ 1). When the spectral indices are kept equal 
(""ad = "iso") the isocurvature component also acquires 
the same spectral index. The multipole dependence of 
the isocurvature contribution to the CMB spectrum can 
now be seen easily from Fig.^c). Isocurvature and cor- 
relation modify more the low multipole end of the spec- 
trum than high multipoles. There is a large difference 
in the power coming from non-adiabatic contributions to 
the Sachs- Wolfe plateau (I < 21) and to the first or sec- 
ond acoustic peak. Hence, the acoustic peak structure 
is distorted significantly leading to a need/possibility to 
adjust it by uj b . However, if the isocurvature spectral 
index is a free parameter it acquires a value that leads 
to a small isocurvature contribution on all scales. This 
happens with will be demonstrated in Fig. ED 

With this large rn so the isocurvature contribution to the 
Ci and to the matter power can remain, e.g., at some 
3.5% level (for our median a — 0.035) on all scales. Then 
the different peak structure of the isocurvature compared 
to the adiabatic one represents only a marginal distortion 
from the pure adiabatic case. This explains why we and 
[l6l ] end up with the "adiabatic value" for ui b . 

The isocurvature effect on the determination of the 
Hubble parameter in our model is more dramatic. With- 
out extra priors we do not obtain an upper bound for h 
(within the analyzed range 0.4 < h < 1.0), see Fig. 01 
and the discussion of the previous subsection. (How- 
ever, with a different choice of the pivot scale we would 
miss the small matter density models and hence obtain 
an upper bound for h. We will discuss this in detail in 
Sec. ED) Applying the SNIa result ft A = 0.70 ± 0.04 to 
get rid of the large h values we get a 95% C.L. region 
0.64 <h< 0.77. 

In Fig. 1111 we compare our model with the pure adi- 
abatic model by showing the C.L. contours in (uj b ,Ho) 
plane. We indicate also the 95% C.L. result u b — 
0.020 ±0.002 from a big bang nucleosynthesis (BBN) cal- 
culation [34( and the Hubble space telescope key project 
result h = 0.72 ± 0.08 0. The 95% region (or even the 



12 



100 

95 
90 
85 
80 
75 
70 
65 
60 
55 



68% C.L. 
95% C.L. 
adi 68% C.L. 
adi 95% C.L. 




We obtain an upper limit (95% C.L.) to the isocurva- 
ture fraction 
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Figure 11: The 68% {solid) and 95% (dashed) C.L. regions 
in the (u)b,Ho) plane for our isocurvature model (black) and 
for the adiabatic model (red). The rectangular gray boxes 
represent a BBN constraint co b = 0.020 ± 0.002 Q and the 
HST constraint H„ = (72 ± 8) km/s/Mpc 0. 



68% region) of our model certainly accommodates the 
HST result, but is only marginally consistent with the 
BBN value of Wf, from [34[ . Actually, the same is true for 
the adiabatic model. On the other hand, concordance 
is achieved with another BBN value uj b = 0.022 ± 0.002 
from [46|. 

For comparison, the similar contours opened up to- 
wards the upper right corner of (u>b,h) plane in [4tH |. 
Moreover, the only "excluded region" was "an upper left 
corner" of their (w&,/i) plane where the BBN and HST 
regions intersected in their Fig. Again we stress that dif- 
ferent data sets were used in [I[| , also neutrino isocurva- 
ture modes were allowed and there was only one spectral 
index. In any case, the considerations in this subsection 
demonstrate that even within "isocurvature models" the 
initial assumptions, e.g. the shape assumed for the pri- 
mordial spectrum, affect considerably the end results. 



D. Isocurvature parameters 

We now turn to the parameters related to isocurvature 
perturbations. The 1-d likelihoods for the 4 independent 
parameters, the isocurvature fraction a, the isocurvature 
spectral index n; so , the adiabatic correlated fraction 7, 
and the spectral index n ac \2 , are shown in Fig. [SI and the 
two derived parameters, the correlation fraction 



a C or = sign(7)vo(I - a) 111 
and the correlation spectral index 

™ad2 + «iso 



lie 



(27) 



(28) 



a < 0.18. 



(29) 



One should be careful about the meaning of this. First, 
a is defined as the isocurvature fraction at our pivot scale 
fco = O.OlMpc -1 . Models with a small isocurvature frac- 
tion at this scale may have a large isocurvature fraction at 
some other scale, depending on how the spectral indices 
for the adiabatic and isocurvature fractions differ from 
each other. Second, since a is defined in terms of the 
primordial curvature and entropy perturbations, it does 
not give directly the relative isocurvature contribution to 
Ci , but that depends also on the shapes of the component 
spectra C*f dl , Cf d2 , C7 ; cor , and Cf°, (i.e., on the transfer 
functions) which depend on the other cosmological pa- 
rameters, and are typically such that the isocurvature 
contribution to the total C; and P(k) is smaller than a. 
Thus the limit to an isocurvature signal in the data is 
actually tighter than would appear from Eq. I|29|) . (See 
Sec. 0) Third, because of the presence of poorly con- 
strained parameters, n a d2 and n- lso in the case of small 
7 or a, the likelihood functions, and thus upper lim- 
its, are sensitive to the priors implied by the choice of 
parametrization. We discuss this last point in Sec. IVII 
Similar caveats apply to the other isocurvature related 
parameters. 

For the "uncorrelated" subset, 7 = 0.0 ± 0.02 the for- 
mal upper limit is larger 



a < 0.22. 



(30) 



in Fig. EJ 



The limit for correlated models is tighter, since the corre- 
lation contribution to the data tends to be larger than the 
isocurvature contribution, due to the transfer functions 
(see Figs. ^ an d HJ), and since |a cor | = y/ a(l — a)\j\ > a 
for small a and moderate 7. 

The isocurvature spectral index has a fairly wide dis- 
tribution covering the range < n iso < 4. The median 
value is ni so = 2.252. The distribution is skew, so that 
the largest marginalized likelihood is at somewhat larger 
values, ri; so ~ 3.0. This peak at ni so comes from the 
large-fiA-models discussed in Sec. IIVBI Otherwise, val- 
ues 1 < rtjso < 3 are preferred. Fig. ^] shows the 2-d 
likelihood for a and rii S0 . 

There are basically two reasons why the data selects 
this range for 7ii S0 . Disregarding the peak structure in the 
Ci spectrum, the overall distribution of power in the data 
over different scales is such that for the adiabatic models 
it favors a scale-independent n a d ~ 1 primordial spec- 
trum. For the isocurvature modes the Ci transfer func- 
tion falls more steeply with k (see Figs. ^ and [21) Thus, 
for the isocurvature contribution not to disturb this over- 
all distribution of power, it needs a larger spectral index. 
The other reason is in the more detailed shape of the 
data. The CMB data clearly does not like the isocurva- 
ture contribution, since it has the wrong peak structure. 
Too small (large) rii so would cause it to show up for small 
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Figure 12: 2-d marginalized likelihood for a and ni s 



Figure 14: 2-d marginalized likelihood for a cor and a. 
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Figure 13: Like Fig.|3c), but now with riiso = 2.252. 



(large) I, even for small a. With k in the middle of the 
data sets, this keeps < rii SO < 4. 

In Fig. El we show the unit-amplitude component Ci 
spectra, for n a d still at 1, but ni so at the median value, 
rii SO = 2.252. Now the effective slope of the adiabatic 
and isocurvature contributions is roughly the same, so 
that the isocurvature contribution is kept low everywhere 
with moderately small a. 



E. Correlation Parameters 

Zero correlation, 7 ss 0, is favored over any other par- 
ticular value for the correlation. However, 61% of the 
models have I7I > 0.1 (Fig. 0. Positive correlations are 
favored over negative ones. 

A strong correlation 7 between the adiabatic and 
isocurvature perturbations however has little effect on 
the data, if the isocurvature perturbations, with which 
the adiabatic perturbations are correlated, are negligi- 
bly small. The signature in the data is better measured 
by the derived parameter a cor (which is restricted be- 
tween ±y/a{l — a) by definition). We see that the 1-d 
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Figure 15: 2-d marginalized likelihood for a cor and n c 



likelihood of a cor is skew (Fig.|SJ); the preference for pos- 
itive correlations that we saw in 7 appears here as a long 
tail towards large a cm . If we add the Gaussian prior 
f?A = 0.70 ± 0.04 to represent SNIa constraints, this tail 
goes away, and the Id likelihood becomes rather symmet- 
ric. Thus the preference for positive correlations is due to 
the large-f^A models discussed in Sec. II V Bl The dip at 
acor = in Figs. El and [TBI does not indicate that un- 
corrected models would be unfavored by the data; rather 
it comes because flat priors for a and 7 lead to a prior 
for a cor which is small for small a CO r- Fig. 1141 shows the 
2-d likelihood of a and a CO r- 

In Fig. El we show the 2-d likelihood of a cor and 
n cm . It shows that positive correlations are connected 
with larger spectral indices n cor than negative corre- 
lations. This feature remains also after applying the 
f?A = 0.70 ± 0.04 prior, although the largest a cor val- 
ues are cut off. For large I7I we have n a d2 ~ 1) an d thus 



(31) 



The reason negative correlations are favored with smaller 
n COI is that then there is a significant isocurvature and 
correlation contribution to the Sachs- Wolfe region of the 
TT spectrum, and the negative correlation now subtracts 
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Figure 16: 2-d marginalized likelihood for 7i a di and n a <i2- 



from it, helping to fit the lowest I WMAP data points 
which lie below the adiabatic spectra (see Fig. I17|) . 

For larger n cor the correlation contribution is insignif- 
icant in the SW region, but becomes important in the 
region of acoustic peaks and for the matter power spec- 
trum. Positive correlations are now favored for the rea- 
sons discussed in Sec. IIVBI This effect remains after 
adding the Gaussian prior J7a = 0.70 ± 0.04, in part 
since this SNIa result favors somewhat larger Qa than 
the CMB+SDSS data applied to adiabatic models. 

Whenever one adiabatic component has negligible am- 
plitude, the corresponding spectral index (i.e., n a d2 for 
7 pa 0, and n a di for |-y | ~ 1) becomes unconstrained (see 
Fig. 116(1 . otherwise it is tightly constrained to be near 1. 
When both components are significant, there is a small 
anticorrelation between n^i and n a( j2 Q, as a red tilt in 
one of them can compensate for a blue tilt in the other 
one, making the sum closer to scale-invariant (as pre- 
ferred by the data). This effect however introduces a 
positive dn/d\nk in the combined spectrum (Sec. lIV^j l. 
which the data does not like, especially with the inclu- 
sion of LSS data, and therefore the anticorrelation effect 
is now more limited than in [3|. 

In Fig. El we show the spectra for our best-fit model. 
This is an example of a low-ni SO , negative-correlation 
model, where the correlation contribution subtracts from 
the Sachs- Wolfe region in the C/. This model has u>t, — 
0.0227, u c = 0.129, 6 = 1.043, r = 0.144, b = 0.962, 
ln(10 10 ^ 2 ) = 3.30, n ad i = 0.988, a = 0.00197, n iso = 
0.388, 7 = -0.67, n ad2 = 0.926, and X 2 = 1459.29. 



Models with Very Large Isocurvature Spectral 
Index 
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Figure 17: The CMB and matter power spectra for our best- 
fit model. The correlated adiabatic component ("ad2", ma- 
genta) dominates over the uncorrelated adiabatic component 
("adl", red) and the negative correlation subtracts from the 
Sachs- Wolfe region of the CMB spectrum. 



We set a very wide allowed range for the isocurvature 
spectral index. While most of the good models had n- mo 
in the range to 4, one of our MCMC chains found an 
apparently very good region where 7ij so was between 5 
and 6. In fact the highest likelihood (x 2 = 1459.20) was 
obtained in this region. We did not have enough statistics 
to assess correctly the relative importance of this disjoint 



good-likelihood region in the parameter space. We dis- 
card this region for reasons explained below. Thus we 
have not included this chain in our full analysis. (And 
therefore we do not take our best-fit model from it.) In 
fact, these models are obviously nonsense, and we discuss 
them just as a warning. 

These models necessarily have a very small a; be- 
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Figure 18: The matter power spectrum of a model with 
«-iso=5.69, which has % 2 = 1459.20. The inset shows the 
SDSS window function for the 16th, i.e. second from right, 
data point. The red markers (*) indicate the values (the- 
oretical matter power spectrum convolved with the window 
function) to be compared to the data points, showing that the 
fit (x 2 ) obtained this way is very good, although the power 
spectrum seems to lie below the data. The problem with this 
kind of matter power is that the window function picks most 
of the "power" from the non-linear regime. Moreover, we 
have calculated the matter power up to k/h ~ 3Mpc _1 only. 
Had we calculated further, the rapidly rising power at large 
k would have caused the red markers (*) to move up leading 
to worse x 2 f° r this model. 



region (Sec lVIfl . 



V. NON-ADIABATIC CONTRIBUTION TO 
THE OBSERVED SPECTRA 

So far we have constrained the non-adiabatic contribu- 
tion to the primordial spectrum in terms of a and a cor (or 
7). Although the isocurvature component can be as large 
as 18% of the total primordial power at our pivot scale 
ko, its role in the observed C; (or matter power) spectrum 
is less significant. This comes because of different behav- 
ior of adiabatic and isocurvature transfer functions as 
discussed in Sees. ITT1 and IIV Dl Moreover, the non-scale- 
invariant spectral index complicates drawing conclusions 
for the observed Cf° and C ; cor from a and a cm , respec- 
tively. Thus we devote this section to finding limits for 
non-adiabatic contributions to the observed spectra. 

We define a relative non-adiabatic contribution to C Z TT 
by 



a 1 



£<TTiso _|_ £<TTcor 



(JTT 



(32) 



where Cf T = C™ 1 



£iTTad2 



£iTTiso 



£tTTcor 



When 

creating MCMC chains we saved this quantity for 1 = 2, 
140, 200, and 700 for each accepted step. By similar man- 
ner we define a non-adiabatic contribution to the matter 
power at k = fej 



cause of the large nj so the isocurvature contribution is 
steeply rising, and only becomes noticeable at the small- 
est scales of our data set. At scales smaller than included 
in our data set, the isocurvature contribution then be- 
comes dominant, and P(k) rises rapidly. 

Thus for most of the data set, these models are es- 
sentially equal to the adiabatic model. The improvement 
over the adiabatic model is then in the "other CMB" and 
"SDSS" data which cover the smallest scales. The fit to 
the SDSS data is however obtained in a rather unnatural 
way. Because the SDSS window functions, that describe 
how the data points relate to the underlying power spec- 
trum, extend to much smaller scales (larger k) than the 
nominal k values of the data points, for these models 
they pick up most of the contribution at these very small 
scales (see Fig. II 81) . Since the perturbations are non- 
linear at these scales, our use of a linear power spectrum 
does not give correct results. (We also suspect that the 
SDSS window functions were not really meant to be used 
for this kind of spectra.) Anyway, these models would be 
ruled out if some smaller scale constraints were added. 

Because our pivot scale is far enough to the left from 
the right (small-scale) end of our data set, these models 
are forced to have a rather small a, which makes the 
measure of this region of parameter space rather small. If 
a smaller pivot scale (larger ko) is used, it becomes more 
likely for the MCMC chains to end up in this questionable 



P iso (h) + P COT (h) 



(33) 



We saved this around the first SDSS data point k\jh = 
0.0154 Mpc -1 and at the last data point kyj/h = 
0.154 Mpc- 1 . 

The range of possible values for ai and a m i is [—00, 1]. 
For example, a; gets negative values whenever Cf OT < 
-Cf °. In the extreme case that Cf = Cf d2 and Cf m = 
— 2C, lso the denominator approaches zero in the absence 



of Cf dl . 



obtained with Cf OT 



On the other hand, the maximum value 1 is 

0. 



Recall that the C; TT s are related to the variance of the 
CMB temperature perturbation by 



1=0 



l)C\ 



TT 



(34) 



We have calculated the Ci for I = 2-2100. In all well- 
fitting models the power at I = 2100 is negligible due to 
diffusion damping. These considerations lead us to one 
more measure of the non-adiabatic contribution 
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Figure 19: Marginalized likelihoods of non-adiabatic contri- 
butions to the observed spectra in our 11 parameter model. 
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Figure 20: Marginalized likelihoods of non-adiabatic contri- 
butions to the observed spectra in "uncorrelated models" 
(7 = 0.0 ±0.02). 



Correlated models. In Fig. 1191 we plot the marginal- 
ized 1-d likelihoods for 012, 0:2001 ami, &mn, and ot- At 
the quadrupole (I = 2) a long tail of a 2 towards nega- 
tive non-adiabatic contribution appears, since the mea- 
sured quadrupole is rather low compared to typical pure 
adiabatic models. The 95% C.L. region spans an inter- 
val —0.46 < CV2 < 0.10. Around the first acoustic peak 
the non-adiabatic contribution is much more constrained, 
-0.024 < a 20 o < 0.079, and at I = 700 the allowed con- 
tribution becomes even smaller, —0.011 < 0:700 < 0.026. 
In the matter power, the limits are —0.06 < a m x < 0.14 
and —0.11 < a m u < 0.51. The latter, quite large values 
come because we do not use any data from larger fc. So 
the spectrum is practically unconstrained after kn. (Re- 
call also our warning example in Figure [TBI ) The likeli- 
hood for the total non-adiabatic temperature perturba- 
tion is quite symmetric with the median at ot = 0.009 
and a 95% C.L. interval -0.075 < a T < 0.075. Hence, 
we conclude that the non-adiabatic contribution to the 
observed temperature perturbation is less than 7.5%. 

Uncorrelated models. Four years ago we [i| found up- 
per limits for an uncorrelated CDM isocurvature contri- 
bution using the first data releases of Boomerang 
and Maxima [H together with COBE data HI. The 
95% C.L. limits were 02 < 0.56 (called a in 6J) and 
"200 < 0.13. Let us update these numbers to reflect the 



dramatically increased accuracy of the data. We approx- 
imate uncorrelated models by applying a Gaussian prior 
7 = 0.00 ± 0.02 when analyzing the chains. Since the 
data does not favor correlation (see Figure [SJ), the sam- 
pling of models with small I7I is very good. For uncorre- 
lated models the correlation component is missing from 
definitions 1(32(1 . (J2HJ, and 135(1 . Then the range for cc;, 
a m i and ar is [0, 1]. 1-d likelihoods are given in Fig. 1201 
The 95% C.L. limits are a 2 < 0.085 and o 2 oo < 0.023. 
So, the allowed isocurvature contribution in the uncor- 
related case has dropped to about one sixth part of the 
limits obtained four years ago. Finally, the allowed total 
non-adiabatic contribution (ot) to the observed temper- 
ature perturbation signal becomes less than 4.3%. 



VI. EFFECT OF CHOICE OF PIVOT SCALE 

When the modes have different spectral indices, the 
relative amplitude parameters a and 7 become dependent 
on the choice of pivot scale kg. In the literature, different 
pivot scales have been used, e.g., fco = 0.002 Mpc^ 1 and 
fco = 0.05 Mpc" 1 , whereas we have chosen an intermedi- 
ate value fco = 0.01 Mpc -1 . 

One can convert the results obtained using one pivot 
scale fco to what one would get with another pivot scale 
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fco, by using the parameter transformation 



a 



(1 - a)(l - ^Dfc^di- 1 + (1 - a)|7|fc n »d2-i + a fcn if 



(36) 



(1 - |7|)fe«adl-l + | 7 |fc«ad2-l 



(37) 



(1 - a)(l - l-fDk"^ 1 - 1 + (1 - a)^^™^ 2 " 1 + afc"' 30 " 1 



(38) 



where A: = fco /fco, and weighting the likelihoods with the Jacobian determinant of this parameter transformation 

£.«adl -l£,n a d2-l£.»Mso-l 



J = 



(1 - |7|)fc"adi-i + | 7 ||;«ad2-i (i - q,)(i - | 7 |)^n a di-i + (1 - a)| 7 |A: n -d2-i + a k n ^~ 



(39) 



This weighting gives the effect of changing from flat priors 
for a, 7, and In A to flat priors for a, 7 and In A. 
Typically we have n a di ~ n a d2 ~ 1, so that 

a a , — - afc™-" 1 (40) 

1 — a + ak niso 1 

7 a 7 (41) 
A 2 « A 2 (l - a + afc™'-" 1 ) - A 2 (42) 

and 

£.n iso -l 

J « ; -2 - fc— 1 , (43) 

(l - a + ak 71 '^- 1 ) 

where the "~" are for small a. Thus, if kg > ho, the 
likelihood of models with large rii so is increased and that 
of small n; so is decreased and the opposite holds if kg > 
k . 

To study the effect of varying k , we both (1) applied 
the above transformation to our results from our main 
run and (2) did shorter MCMC runs (8 chains) using 
fc = 0.002 Mpc -1 and k = 0.05 Mpc" 1 . Both methods 
should give the same result if the MCMC runs have suffi- 
cient statistics. In practice, the results were close to each 
other for k = 0.002 Mpc -1 , but for fc = 0.05 Mpc -1 , 
our original run had insufficient sampling at large rii so , 
for the reparameterization to give meaningful results. 
We show in Fig. [2] the resulting marginalized likeli- 
hoods for the (primary) parameters most affected. For 
fco = 0.002 Mpc" 1 the result shown is by method (1), 
but for fco = 0.05 Mpc" 1 by method (2), since it had 
better statistics. However, these results should only be 
taken as indicative, especially for fco = 0.05 Mpc" 1 as 



the statistics was not nearly as good as in our main case, 
fc = 0.01 Mpc" 1 . 

The 1-d likelihoods of u>b, r, b, and n^i did not change 
significantly. Thus these parameters are not sensitive to 
the choice of pivot scale. The parameter affected the 
most is 7ii so , where we see very clearly the shift to smaller 
(larger) n,i SO as the pivot scale is increased (decreased). 

Consider first the change to a large pivot scale fco = 
0.002 Mpc" 1 , corresponding to Iq ~ 28. Now a lot of 
weight is given to models with a "red" isocurvature spec- 
trum, riiso < 1- For these models the isocurvature con- 
tribution is significant in the SW region of the CMB 
spectrum, and negligible elsewhere. Accordingly, neg- 
ative correlation 7 is favored, since it subtracts power 
in the SW region where the data is below the adiabatic 
model prediction. A red correlated adiabatic index n a d2 
is favored, as the \ow-l boost in the negative Cf or tends 
to win over that in the positive Cf A2 . Because of the 
negative correlation contribution, somewhat larger am- 
plitudes A are favored (not shown in Fig. l21fl . With very 
little weight at ri; so ~ 3, the large-l^A models are elimi- 
nated, so the tails in the to c and 9 distributions disappear. 

The 1-d likelihood for the isocurvature fraction a is 
surprisingly close to the fco = O.OlMpc" 1 case. Thus our 
upper limit a < 0.18 seems to be more robust than one 
might have thought, and applies over a fairly large range 
of scales. 

Consider then the change to a small pivot scale fco = 
0.05 Mpc -1 , corresponding to Iq ~ 700. This has the ef- 
fect that the problematic "high likelihood" region around 
n- mo ~ 5-6, discussed in Sec. IIVFI acquires a much 
larger measure in the parameter space, increasing the 
marginalized likelihood of these n- lso values. These mod- 
els have now a large weight in the 1-d likelihoods of 
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Figure 21: 1-dimensional likelihoods for w c , a, n[ so , 7, 
and n a d2, using three different pivot scales, fco = 0.002Mpc _1 
(dashed, red), fco = O.OlMpc" 1 (solid, black), and fco — 
0.05MPC" 1 (dotted, blue). 



all parameters. While they had a very small a (fco = 
O.OlMpc" 1 ), they have a rather large a(k — 0.05Mpc~ ) 
(see Eq. JH5J), and thus the a distribution now extends 
to large values. At the 95% C.L. we obtain a(k = 
0.05MPC" 1 ) < 0.56. The "bump" in the uj c distribu- 
tion around ~ 0.11 is also due to these models. The way 
the SDSS window functions collect power from smaller 

" O fpisor 



scales fSec. [IV Fp of P lso (fc) allows a smaller "shape pa- 
rameter" VL m h for P ad (fc), which leads to the smaller u> c - 
As explained in Sec. IIVFI we do not take these models 
with n; so ~ 5-6 seriously. Thus the fco = 0.05 Mpc" 1 
case here should just be taken as a warning for what 
may happen with extreme values of spectral indices in 
this kind of studies. 

In general, the pivot scale should be chosen to be in 
the middle of the data set. If fco is near either end of 
the range of scales covered by the data, the spectral in- 
dices of components which are subdominant in the main 
data region become unconstrained in the direction which 
causes this component to blow up outside, or near the 
edge, of the data set. 



VII. COMPARISON TO BELTRAN ET AL 

Due to the many differences in approach discussed in 
Sec.|n]tlie comparison of our results to those of is not 
straightforward. If we include an Oa = 0.70 ± 0.04 prior 
to mimic their use of SNIa data, the (large 9, small u> c ) 
models we found but they did not, disappear from our 
results. The remaining minor differences in the determi- 
nation of standard (adiabatic) cosmological parameters 
are mainly due to a different choice of pivot scale fco- 
When we shift to the same pivot scale, fco = 0.05 Mpc" 1 
they used, our results approach theirs. Our results for 
this pivot scale are however contaminated by problem- 
atic very large n; so ~ 5-6 models, whereas they have 
imposed an upper limit nj so < 3, so the results are not 
directly comparable even in this case. 

The parameters related to isocurvature perturbations 
are defined with respect to the pivot scale. Our upper 
limit to the isocurvature fraction a at pivot scales fco = 
0.002 Mpc" 1 and fc = 0.01 Mpc" 1 is much tighter than 
their upper limit of about 60% at fco = 0.05 Mpc" 1 . Our 
unreliable fco = 0.05 Mpc" 1 upper limit a < 0.56 agrees 
with that limit. 

Because of different choice of correlation parameters, 
the results for correlation are best compared in terms 
of our o w, w hich equals their —fly/ a(l — a) plotted in 
Fig. 3 of jig . Our discovery of a preference for positive 
correlations at large n; so is in agreement with their result 
(with the opposite sign convention). 



VIII. DISCUSSION 

We have used CMB and large-scale structure data 
to constrain models where the primordial perturbations 
have both an isocurvature and an adiabatic component, 
allowing for different spectral indices for these compo- 
nents, and a possible correlation between them. We re- 
stricted these models to a spatially flat (f2 = Q\ + Q, m — 
1) background universe. 

The basic conclusion is that the data clearly disfavors 
the presence of isocurvature perturbations. This makes a 
likelihood study of such models problematic, since once 
the isocurvature contribution is small, the related spec- 
tral indices become unconstrained. When some of the in- 
dependent parameters are unconstrained, the likelihood 
function becomes sensitive to the implied prior due to the 
parametrization used. We demonstrated this by chang- 
ing the pivot scale used to define our isocurvature and 
correlation fraction parameters. 

The problem with spectral indices does not occur when 
a model has only one independent spectral index. It 
would also not occur if the data would clearly favor a 
nonzero fraction for any component whose spectral in- 
dex we have as an independent parameter. 

Perhaps a better parametrization of isocurvature mod- 
els would be to use the amplitudes at two different scales 
(e.g. at fc m i n and fc max ) as the independent parameters 
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for the likelihood analysis, instead of an amplitude at one 
scale and a spectral index. The spectral index would then 
become a derived parameter. We suggest that one tries 
this approach in future studies, since it might: 1) Lead 
to a much faster convergence of the MCMC chains be- 
cause the unconstrained spectral indices would be miss- 
ing. 2) Remove a possible bias towards zero isocurvature 
amplitude models, which was a result of blowing up the 
parameter space volume upon marginalization caused by 
unconstrained n; so in case of small a. 3) Prevent the 
feature that with too large kg the integration measure 
(weight) of models with extremely large nj so becomes ar- 
bitrary large. 

For models with the largest isocurvature fractions at 
the pivot scale ko = O.OlMpc" 1 , which is roughly in the 
middle of the data set used, the isocurvature spectral in- 
dex is constrained to be in the range 0.5 < n- lso < 3.5 
which prevents the isocurvature contribution from ris- 
ing too high either in the small- or large-scale ends of 
the data used. If one moves the pivot scale to smaller 
(larger) scales the upper (lower) limit to n iso is relaxed, 
or removed, as the rising part of the isocurvature spec- 
trum moves outside the data range. 

Of the standard (adiabatic model) cosmological pa- 
rameters, the determination of the baryon density LOb, 
the primordial perturbation amplitude A, the adiabatic 
spectral index n a( j, the optical depth due to reionization 
t, or the bias parameter 6, is not significantly affected by 
a possible isocurvature contribution. On the other hand, 
models with a smaller CDM density lo c and a larger sound 
horizon angle 9 become acceptable. This means that we 
cannot even rule out models with Hq > 100 km/s/Mpc 
and n m < 0.1 (at 95% C.L.) using CMB and LSS data 
alone. 

We obtained an upper limit a < 0.18 (95% C.L.) for 
the CDM isocurvature fraction for models where corre- 
lation is allowed between the isocurvature and adiabatic 
contributions. This limit is somewhat tighter than the 
corresponding limit for uncorrelated models, since cor- 
relation causes a stronger signature in the data than an 
uncorrelated isocurvature perturbation. 

Here a is defined as the ratio Vs/^Pn + Vs) 01 the pri- 
mordial entropy and curvature perturbation power spec- 
tra, at a pivot scale fcoj an d our upper limit applies 
for both k = 0.002 Mpc" 1 and k = 0.01 Mpc" 1 , and 



presumably also for the range in between. The value 
a = 0.18 corresponds to / lso = {V s /Vn) 1/2 = 0.47. For 
smaller scales (larger k) our results are less conclusive, 
since there the constraint on a relies more on the large- 
scale structure (SDSS) data, whose use is problematic 
for a steeply rising (large rii so ) isocurvature contribution. 
However, our results for fco = 0.05 Mpc -1 are not in dis- 
agreement with the upper limit of 60 % for a obtained 
in |16( using this pivot scale. 

In the observed temperature anisotropy signal the 
amount of non-adiabatic contribution is \oit\ < 0.075 
at 95% C.L. in our correlated isocurvature model. The 
upper limit becomes tighter in the uncorrelated case, 
a T < 0.043 at 95% C.L. 

In models with a large isocurvature spectral index, 
rij SO ~ 2-4, a positive correlation between the adiabatic 
and isocurvature perturbations is favored. The correla- 
tion contribution appears then in the acoustic peak re- 
gion, where the effect of a positive correlation is to shift 
the acoustic peaks towards larger multipoles /, which 
then favors a larger sound horizon angle 9 to push the 
peaks back to where the data has them. To satisfy also 
the large scale structure data, smaller CDM densities ui c 
are then favored. These effects translate into a larger Hq 
and a smaller Q m (larger 

In models with a small isocurvature spectral index, 
"iso ~ 0-2, a negative correlation is favored. Here 
the correlation contribution appears in the Sachs- Wolfe 
region, where this negative correlation brings the Ci 
down to better agree with the small large-scale CMB 
anisotropy seen by WMAP. 
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